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ABSTRACT 

A laminate composite containing an orthotropic layer with a 
crack situated normal to the interfaces, and bonded to two ortho- 
tropic half-planes of dissimilar materials, is considered. The 
solutions for two different classes of orthotropic materials are 
presented separately. In each case, the problem is first reduced 
to a system of dual integral equations, then to a singular integral 
equation which is subsequently solved numerically for the stress 
intensity factors at the tip of the crack. The effect of the 
material properties on the stress intensity factor is investigated. 
Two cases, the generalized plane stress and the plane strain, are 
treated simultaneously. 
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1 . 


INTRODUCTION 


The increased use of composites can probably be attributed to 
the fact that they are a group of structural elements which opti- 
mize the requirements of safety and economy. As a result, this 
has prompted the solution to many problems involving various com- 
posite geometries [l]-[4], [10], [13] -[17], One particular 
composite geometry of practical importance is the multi-layered 
isotropic medium consisting of many layers, including the presence 
of possible flaws. This problem has been treated by Hilton and 
Sih [1] by simplifying the geometry to the form of a single layer, 
with an internal crack normal to the interfaces, bonded to two 
dissimilar half-planes with averaged elastic properties. Through 
the application of an integral transform technique, a Fredholm 
integral equation is obtained and solved to obtain the crack tip 
stress field. Recently, the same problem has also been considered 
by Bogy [4] . Furthermore, the case of a crack touching the inter- 
faces has been investigated by Ashbaugh [2] and Gupta [3] , where 
the stress state around the crack tip has been found by the method 
of integral transforms, the singularity having been shown to be 
different than 1/2. 

In this paper, the solution of the problem treated by Hilton 
and Sih will be extended to the case of orthotropic materials. 
Since the orthotropic materials can be classified in two groups, 
there are four possible material combinations. Two cases, where 
the layer and the matrix are both of the same material type, are 
formulated, and the corresponding solutions are presented. The 
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solutions to the other two combinations can be obtained, in an 
identical manner. In each case, the problem is formulated in' 
terms of integral transforms, using the procedure described by 
Sneddon [5] . The dual integral equations resulting from the mixed 
boundary conditions are solved by reducing them to a singular 
integral equation, and then employing the numerical technique 
described in [7]. For certain choices of material parameters, 
the effect of crack size on the stress intensity factor is shown 
in Figures 2 - 6 . 
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2 . 


GENERAL FORMULAS 


For an anisotropic elastic body the equilibrium equations, the 
stress-strain and the strain-displacement relations in the absence 
of body forces are: 
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Substitution from (2.3) into (2.2) and (2.1) yields the equilibrium 
equations in terms of the displacements: 
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(2.5) 


Defining the inverse of the matrix by 

[a] = [A]“^ 

and noting that for orthotropic materials 
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which reduce to the Navier's equations for isotropic materials. 
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Since the problem under consideration is one of the plane 
theory of elasticity, the formulas of the plane strain and the 
generalized plane stress will be given below. 

2.1 Plane Strain 


From the conditions of plane strain, i.e., 
w = 0 , u = u(x,y) , V = v(x,y) 
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Furthermore from (2.2) we conclude that for the bodies possessing 
at least one plane of elastic symmetry (the plane perpendicular to 
z-axis) 
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and the first two equations of (2.4) can be solved to give the 
displacements u,v, the third being automatically satisfied. 


In particular, for orthotropic bodies we have 
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2.2 Generalized Plane Stress 
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Substitution from (2.3), (2,13) and (2.14) into (2.1) yields 
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For orthotropic bodies we have 


[a] = 


1 

V 

Z2£ 

E 

y 

0 


1 

^yx 

0 

E 

x 


E *A 

y 

Ey-A 

V 

_ ^ 

E 

X 

1 

0 

, [A] = 

V 

xy 

1 

0 


E • A 
x 

E ’A 
X 

0 

_ 

0 

l/G 

xy_ 


0 

0 

G 

xyj 


A = 


E *E„ 
X y 


(1 - V V ) 
' xy yx 


and (2,16) reduces to 


A ^ + A 
11 2 33 


92u 

9y2 


+ (A^o +A-i<x) 


2 

9^v 


33 ^"12' 9x9y 


= 0 


^ V + A LJL + (A +A — = 0 
'33 ^..2 ^ ^22 ^..2 ^ '^33 ^ 12 ^ 3x9y 


(2.17) 


(2.18) 


9x“ ““ 9y 

Hence for the orthotropic materials equations (2.11) and (2,18) can 
be expressed as 
.2 


3 


1.2 

9x 


1 ^. e 

3x2 


92u 

®3 

92v 

9y2 

3x9y 

9y 


9 2u 
9x9y 


= 0 


0 


(2.19) 


9 



where 


A 
^ A 


11 


66 



At ^ 

g 1 + for plane strain 

^ ^66 




1 + 



for generalized 
plane stress 

( 2 . 20 ) 


3. formulation of the problem 

A single layer of width 2h is assumed to be bonded to two half- 
planes of dissimilar materials. Both the layer and the matrix are 
considered orthotropic. A crack of length 2a is symmetrically 
located being normal to the interfaces (see Figure 1) . The case 
a<h will be studied. The case a=h is given in [10]. The crack 
surface loading consists of normal stresses distributed symmetri- 
cally. The solution of the skew- symmetric case can be obtained 
similarly. A superscript * for the elastic constants will refer to 
the matrix. 

The approach described in [5] will be used. First, we will 
obtain the solutions of (2.19) satisfying certain boundary condi- 
tions of the layer and the matrix. The combination of these solu- 
tions will be forced to satisfy the rest of the boundary conditions. 

3.1 Solutions u^(x,y) , v^(x,y) 

Assume 

oo 

^q(x,y) = f J f(t,x) costydt 
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(3.1) 
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Hence u^(x,y) and v^(x,y) are determined in terms of unknown func- 
tions A(t) , B(t), C(t), D(t). Note that 3’s here correspond to the 
layer. 
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3.2 Solutions U 2 (x,y) , V 2 (x,y) 


Assume 


U 2 (x,y) = - ! f^(t,x) cos tydt 
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Solving (3.8) we obtain 
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and the other constants are as in u^(x,y) and v^(x,y) 


3, 4 Stress and the displacement fields of the laminate composite 

The displacements of the layer and the matrix will be expressed 
as 

u^(x,y) = u^(x,y) + U 3 (x,y) 

v^(x,y) = v^(x,y) + (x,y) for the layer, and 

U 2 (x,y), V 2 (x,y) for the matrix (3.11) 


From the symmetry conditions we have 
B(t) = -A(t) , D(t) = -C(t) 

Furthermore from the condition that U 3 (x,y) and V 3 (x,y) will 
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tend to zero as the unknowns G(t) and F(t) are eliminated. 

Also, the condition that U 2 (x,y) and V 2 (x,y) go to zero as 
will be used to eliminate K(t), M(t) . 


The stress field is obtained from (2.14) and (2.17) for the 
case of generalized plane stress. Stress field for the plane 
strain may simply be obtained by substituting ^j^2"^^ll' ■^12'^‘^22' 
1/A^^ and 1/^22 place of ^xy' 

stress expressions. 

In Appendix A, it is shown that s^ and S 2 are either real or 
complex conjugates. Accordingly, the expressions below are given 
for two different types of orthotropic materials. Hence, after 
intermediate manipulations we have. 
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For Material Type II 
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where A(t) , C(t), E(t), H(t) , K(t), M(t) stand for the combinations 
of the previous unknowns. 

The bielastic constants and are given in Appendix B. 

Hence, there are four possible material combinations for the 
laminate composite. However, we will obtain the solutions for two 
cases where the matrix and the layer both are either of material 
type I or II. The other two combinations can be derived in an 
identical manner. 

The boundary conditions of the problem may be summarized as 
follows . 

At X = h 

Uj^(h,y) = U 2 (h,y) ; v^(h,y) = V 2 (h,y) 

C35i(h,y) = ^ '^xy2^'^'y’ (3.16) 

At y = 0 

'^xyl^^'^^ = 0 , |x| <h; 

V 2 (x,0) = 0 , 1 x| > h 

cf^^(x,0) = -p(x) , |x| < a; 

At X = 0 

u^(0,y) = “^xyl^^'^^ " ® 


= 0 / |x| > h 

(3.17) 

v^ (x,0) =0 , a<|x|<h (3.18) 

(3.19) 
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4. 


THE DERIVATIOIsr OF THE INTEGRAL EQUATION AHD THE SOLUTION 


4, 1 For Material Type I 

The stress and the displacement fields given by (3.12) and 
(3.13) automatically satisfy the boundary conditions (3.17) if we 
choose 

H(t) = - E(t) (^-1) 

^12 


On the other hand the mixed conditions (3.18) give 

TT (1 - V V _) 

^ xy yx^ (X 0) = 

2E„ V' ' ' 


» -|u),|ty//67 Yii “Iwolty/ZST 

lim / E(t) (Y.e ^ ^ - Yg e ) cos tx tdt 

y-^0 o ^ ‘ ® ^12 


+ / [2y^A(t) cosh(o)^tx) + 2ygC (t) cosh (w^tx) ] tdt 

I x| < a 


ir (1 - V V ) 

xiLJ2L_p(x) , 


2E 


Vj^(x,0) = - I Yi 9 / E(t)cos txdt = 0 , 


a < X < h 


(4.2) 


To obtain the proper singularity, we will follow the procedure 
described in [13] and [14], and define a new unknown function by 


3v. (x,0) 2 7 , 

— i- — = — Ynn J E(t)sin tx tdt = c|)(x) 

o 

such that 

<|)(x) = 0 , |x| >a 

and finally invert (4.3) to get 


(4.3) 


Y_-.tE(t) = J (j) (x) sin tx dx 
o 


(4.4) 
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Substituting from (4.4) into the first equation of (4.2) one 
arrives at the singular integral equation. This will be done by 
making use of the integrals given in Appendix C. Hence/ 
a °° 

Y 2 q / / [2y^A(t) cosh(w^tx) + 2ygC (t) cosh (w^tx) ] tdt 

-a o 


IT (1 - V V ) 

' 


< a 


(4.5) 


where the definition of 4>(t) has been extended into the (-a,0) 
range as an odd function. The single- valuedness of the displace- 
ments can then be expressed as 
a 

J <l> (t) dt == 0 (4,6) 

-a 


The next step is to determine A(t) and C(t) from the equations 
(3.16). From (3.16), (3,12) and (3.13) after some manipulations 
and an inversion we obtain 

-2A (t) sinh (oj^th) - 2C (t) sinh (w^th) 

- I I th - I ti)^ I th 

+ L(t)e ^ + N(t)e 

= — J cos tydt / E(C) [o 
o o 


^11 -\^n\^Y/^ 


Y 


] sin d^ 


12 


= I / E(C)sin ?hdC [— 2"^ ^2 “ 2 

o (0)^C /B5) + t^ 


I ^1^3 U/ 

^ ^12 (w^C^/35) + t' 


] 


= Rj^(t) 
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2 A{t) 3 ^cosh(Wj^th) + 2 C (t) 3 gCosh 


, * -|io*|th * * -hjlth 

+ sign (w^) 3 ^L(t)e + sign (oj^) 3 gN (t) e 


= — / sin ty dy / E(^) [sign(w^)3ge 
o o 

Y 




11 -kolCy/v^ 

— sign (w^) ] cos ChdC 


12 


00 ^ 

= ^ / E(C)cos 5 hd 5 tsign(w, )Bq — 5—5 ~ 

n ^ ^ (oifr/e^) + 

*- 1 = R.(t) 


r 


11 

12 


sign(w3)eio -— 2 ^ 


(M35 /S5) + 


* 

’l' 


- 1 w, I th * - 1 w- th 


[-2Y3A(t) cosh(w^th) - 2Y4C (t) cosh (w^th) 
* 

+ Yi7Yi3L(t)e 

00 

- ^ f cos tydy / E(^) [Y^^e 


+ Yi7Yi4N(t)e 

- Ujl ^y/’^ 


] t 


Y 


- Y 


2 Y 


11 -^3lCy/v^ 

— e 


] cos Ch Cd^ 


12 


9 1 1 ^^5 

^ f E(C)cos ^h^d^[Yj^ — 


733) + t' 


“ Y 


11 


hoU /717 


^ ■^12 ( 0 )^^^/$^) + t 


J ] = tR^(t) 


4 
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[-2TQA(t)sinh((0j^th) - 2 Yj^qC (t) sinh (lo^th) 




* 

- 0)^ th 


= ^ / sin tydy / E(C)Y;L 3 ^[e 


- 1 I 5y//B^ 


- 10 ) 3 ! th j ^ 


-l“,UY//e7 

- e ^ ] sin 


2 Y ^ 

/ E(5)sin ?h 5d?[ —^2 


tr 


= tR4(t) 


(( 0 ^ 5 Veg) + (<03^^65) + 


(4.7) 


Functions R^Ct) can be obtained in terms of (j> (t) by substitut 
ing from (4.4) into (4.7) and using the integrals given in Appendix 
C. Hence ^af ter lengthy manipulations, 

a -(h-n)t/^/|w I 

f [p 

2 Yn^t i ^ 


Rl(t) 


19" -'a 


Yii ~(h-n) t/^/ jw^l 


5/ j-3 


Y 


12 

a 

/ 

19" -a 
Y' 


] (J>(n)dn 

- (h-n ) t/B ^/ 1 


^2<t) = 2^ _l [sign(<*>i)S9e 

- (h-n) t/^/| 


Y 


11 

12 


sign(w2) 


] 4> (n) dn 


R3(t) 


a - (h-ri) t/B^/l ( 0 ^ 

2y77t [ 


19 -a 0 ). 


^11 ^2 
^12 |t0o 


- (h-n) t/B7/! w 


5' ’ 3 


] {j)(n)dn 
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R4<t) = 2y 






19" -a 


w. 


- (h-n) t/3T/| w. 


] (f){n)dn 


(0. 


Solving (4.7) for A(t) and C(t) and substituting into (4 
gives the following singular integral equation. 


I ; + / k(x,t)<(.(t)dt = - 

~a ~a 


“ ;y^ -- p<-> ■ 

'20 y 

lx| < a 


subject to 
a 

J (J) (t) dt = 0 
-a 


where 

«> - (h-t)n/^/|‘^-i I 

k(x,t) = — — J [k. (x,ri)e 
^^20 o ^ 

-(h”t)r|/BT/l(02l 

+ k 2 (X,n) e 1 dri 

kj^ (x,n) = 2 y"W t i T [•V 7 f 5 (x,n)f,(ri) + Ygf g <^' 1 ) f g (n) ] 
k^Cx^n) = 2y^g^f (n) + Ygfg (x,n) f (ti) ] 


and 

f(n) = f 3 (n)f 3 ^(n) - f2(n)f4(n) 
fj^(ri) = Y25 + Y 24 't^^^ 
f 2 (Tl) = Y 29 

f3(n) = Y27 ‘Y 26 tanh (oj^^nh) 


(4.8) 
.5) 

(4.9) 
(4.10) 


(4.11) 
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f4(n) = Y23 Y 24 tanh(o)^nh) 

f^(n,x) = cosh (Wj^nx) /cosh (Wj^nh) 
fg(n,x) = cbsh (co^Tix) /cosh (co^nh) 


f^ (n) 

" ^38 

)■ Yg^tanh ((D^rjh) 

fg (n) 

0 

11 

+ Y 42 ^tanh (w^rih) 

fg (n) 

^42 

+ Y43tanh ((jOghh) 

fio(n) 

= X44 + Y45tanh(w^nh) 


(4.12) 


Changing the variables 

t = aT , X = ax and 

(j>(t) ^ f p(x) = Pq(x) 

k^(Xr''^) = ak(x,t) 

(1 - ^xv^vx^ 

g(x) = - 2 y E ^ ' Po^^^ 

we obtain 

i / dT + J k^(x,T)4>o(T)dT = g(x) , 

^ 1 1 


/ (j)^{T)dT = 0 


(4.13) 


xl < 1 


(4.14) 


The solution with the proper singularities will be sought in 
the form 

^ (T) = MlL- (4.15) 

° /I - 

where F(t) is Holder continuous in the interval 

Hence following the procedure described in [7] we arrive at 
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N F(T.) 

:=1 j i 


-n- 

I iS- F(T^) = 0 


k=l “' 


(4.16) 


Tj = cos^ ( 2 j-l) , 


j — 1 /.../N 


TTX 

X . = cos^ 
N 


X — If. .of N“ 1 


(4.17) 


From (4.16) N unknowns F(t^), j =l,,..fN can be solved. 


4.2 For Material Type II 


Choosing 


H(t) — ~ TT E(t) 

11 


(4.18) 


we see that the expressions (3.14) and (3.15) for the stress and 
the displacement fields automatically satisfy the boundary condi- 
tions (3.17). 

The mixed conditions (3.18) give 
I'd - %v''vx> 

®yld-0) - 

CO 

J { 2 A(t) [-A^sin(w 2 tx) sinh ( 0 )^tx) + A^cos (ca^tx) cosh (w^tx) ] 
o 

+ 2C (t) [A^sin (a> 2 tx) sinh (w^tx) + A^cos (w 2 tx) cosh (w^tx) ] } tdt 

00 “ 2 ^^ ^^ 2 ^^ 

+ lim f E(t) [(A-cos + A^sin ) 

y -^0 o /F 5 /ST 


Woty 


- T— (-A0COS— — + A_sin— — ) j 

*11 ® / 8 i ’ / 6 ? 

ird - V V „) 

P(x) , |x| <a 


2 ty -|o)j^|ty//6^ 


COS txtdt 
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a < |x| < h 


(4.19) 


v^(x,0) = f ^16 / E(t)cos txdt = 0 , 


Now, in a similar way, we will define the new unknown by 

(x,0) 9 7 /V (A on\ 

i = _ ± f E(t)sin tx tdt = 4> (x) (4.20) 

Ox TT 16 ^ 

such that (f>(x) = 0 for |x| > a. 


Inverting (4.20) we obtain 
a 

-A-^tE(t) = / (x) sin tx dx 

o 

Substitution from (4.21) into the first equation of (4.19) 


yields 

cl ^ 

J + / {2A(t) [-A sin(w 2 tx) sinh(Wj^tx) 

^ -a ^ o 

+ A^cos ^ 

+ 2C (t) [A^sin (^ 2 '^^^ A^cos (w 2 tx) cosh (o)^tx) 1 } tdt 

^ _ in xy yx^. ^ 1^1 (4.22) 

2Ey 

where the extension of the definition of i)>(t) into the (-a, 0 ) range 
requires the single-valuedness of the displacements, i.e., 

J (J>(t)dt (4.23) 

-a 

A(t) and C(t) are determined from (3.16). 

From (3.16), (3.14) and (3.15) after some manipulations and 
an inversion we have 


27 



2 A(t) cos (a)2th)sinh(Wj^th) - 2 C (t) sin (0)2 th) cosh (w^^th) 

- 1 I th . * , , - 1 <^3^ 1 th 


+ K(t) cos (W2th) e 


00 

= i ; E(£)sin 5 hd? 


28 ^12 . 

? A t 


+ M{t) sin(u) 2 '^^^ ® 

A 






11 




11 


U 3.5 2 ( 0 ^?^ (*),? 

4— + ( 4_ + t) 4— + ( -= - t) 


= H^(t) 


5 /85 ^5 

-2A(t) [ 3 ^cos ({02th) cosh (ijJj^th) - g^sin(w2th) sinh (w^th) ] 

- 2C(t) [e!^cos<W2th)cosh(w^th) + 3^sin (W2th) sinh (w^^th) ] 


- 10 )^ th 


+ K(t) [-By sign ((0^) cos ((02th) - By sin(w2th)]e 
+ M(t) [By cos(u)2th) - By sxgn{Wj^)sxn(0)2th) ]e 


00 

= i / E( 5 )cos ShdC 


“33 


5 + 


f*. 
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( 0 ? 5 ^ «o 5 2 W? 5 ^ U,€ 2 

4- + (^+t) 4 — + (-^-t)^ 

Pc ./S- Pr /g- 


= « 2 <*^> 


■5 /S^ ►'S 

{- 2 A(t) [-A 2 sin(ti 52 th)sinh(in)j^th) + A^cos (0)21*^) cosh (oij^th) ] 

- 2 C(t) [A,sin((D2th)sinh(c0j^th) + A2Cos(M2th) cosh(oij^th) ] 

* :S- * 

- K(t) [A23^sin(w2th) + A22COS (W2th) ]e 

* * -|“i|th 

+ M(t) [A^^cos (032th) - A^^sin(o)2th) ]e }t 


00 

= i / E( 5 )cos £hCdS 

IT ^ 



A t 

^ 18 ^ 

2^2 
0 ), g 

OJ.? 2 

X 

+ (— + t) 

L ®5 



hBl 


03 ^^ 

^ - t) 


3 . 




= tH^(t) 
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{- 2 A(t) [A^cos(ai2th)sinh((o^th) - A^^sin (co^th) cosh (w^th) ] 

- 2 C{t) EAgSin(oo2th)cosh(u)^th) + A^^cos (w^th) sinh (oj^th) ] 

* * , -|ui*[th 

+ K(t) [A^^cos (oj^th) + A^^sinEw^th) ]e 


+ M(t) [A 24 sin(w 2 th) - A25COS {o)2th) ] e 


it . 


00 

= ^ / E( 5 )sin?hSd? 


19 


“ll5 


2^2 

“l^ 

TT 


+ ( -f t) 



l“j.U 

^19 



2 

X 

+ ( — + t) 

/FT J 



tH^(t) 


( 4 . 24 ) 


From ( 4 . 21 ) and ( 4 . 24 ) and the integrals given in Appendix C 
H. (t) can be obtained in terms of (j>(t), i.e.. 


A a 

Hi(t) = J {h^^cosi (1)2 (h-n)t] 


-a 


- 0). (h-n)t 

+ sin[ (h-n)t]}e (l)(n)dn 


H2(t) = ^ 


f { COS [ ^2 


-a 


H. 


- u) , (h-n) t 

+ A^gSin[ W2 (h~ri) t] }e <(5(ri)dTi 

A^- a -0), (h-n)t 

(t) = / sin[w^ (h-n) t]e (|>(n)dTi. 

-a 


H (t) = ~ cos[ oi^th-n)!] 

^ ^ -a 1 


- 0 ). (h-n) t 

- sin[ a>2 (h-n)t]}e ^ 4 >(n)dn ( 4 . 25 ) 
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Solving ( 4 . 24 ) for A(t) , C(t) and substituting into ( 4 . 22 ) we 
obtain the following singular integral equation 


1 r . f 

ir i t -X •' 

— 3 . —a. 


6 (x, t) (j) (t) dt 


subject to 

'5 

a 

/ 4>(t)dt = 0 
-a 

where 


~ ^xy^yx* 

2A2,Ey PW ' 


X < a 


( 4 . 26 ) 


9 (x,t) = / {ej^(x,n)cos [ u>2 (h-t)n] 

27 o 

- w, (h-t) Ti 

+ 02 (x,n) sin [ 0)2 (h-t) nl }e dr\ 

e^(x,n) = g [g^ ^^^99 ^x,n) + g^ (h)gj^Q(x,n) ] 

02(x,n) = y tg^ g (x,n) + g 3 (n)g^Q(x,n) 3 ( 4 . 27 ) 

g^(n) = gjL(n)g4(n) - g2(n)g3(n) 

A A 

g2(n) = ~ ^ cos (i02nh) + [— sin(w2rih) + cos (W2nh) ] tanh (Wj^nh) 

A A 

g2<n) ^ cos(o)2nb) + [ ^ tanh(o)^nh) + 1] sin ((02rih) 

93(^1) = “ A^^cos (o)2Tih) + [A^^tanh (w^nb) - 1] sin ((j02nh) 

g^Cri) = “ A^^cos(o)2hh) + [-A ^^sin ((02r|h) + cos (t02nh) ] tanh (w^nh) 

( 4 , 28 ) 


gs (n) = A52gi(n) + A33g2(n) + A53g3(Ti) + 

g^U) = A|^gj^(n) + A^5g2(n) + Ag^g3(n) + A^gg^(n) 

gg (Ti) = A.^^gj^(n) + A^^g2(r|) + A^3g3(ri) + A^^g^(n) 

30 


( 4 . 29 ) 



sin (w 2 nx) 


sinh (Wj^nx) 
(x,ri) - cosh (Wj^nh) 


gj^Q(x,n) = 


cosh (w^nx) 


SS' shTM^^ph) <=°s(o.2nx) 

Changing the variables 

t = ax , X = ax 


(j)(t) = 4 >^(t) , 


^Q(XrT) = a0(x,t) 


p(x) = Pq(x) 


(1 - V V ) 

P„(X) 

27 y 

we then have from (4,26) 

, 1 <)> (t) 1 

¥ _[ -rnr + _/ eo(x,T)4>„(T)dT 


/ (f) (x)dT = 0 
-1 ° 


= r(x) , I x| < 1 


(4,30) 


(4.31) 


(4.32) 


Expressing the solution as in (4.15) we arrive at 
N F(t.) , 

, I — N t T. - X- ■^®0<Xi'Tj) 1 = r(Xi) , i = 1 N-1 

k=l 1^1 


N 


I ^F(x ) = 0 

k=l J 

where x. and Xj^ defined by the equations (4.17). 


(4.33) 


5. STRESS INTENSITY FACTOR 

Due to the symmetry x^^j^(x,0) - 0, hence the stress intensity 
factor will be defined by 

K = lim /2 (x-a) ct . (x,0) (5.1) 

x-^a ^ 
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5.1 For the Material of Type I 

From the equation (4.9), the dominant part of ayj^(x,0) which 
contributes to the stress singularity at the crack tip can be 
expressed as 


a , (x,0) = (x,0) 

yl ' 7 t( 1 - V V ) •’ t- X yl ' ' ' 

^ xy yx' -a ^ 


where from (4.15) 


(5.2) 


0 (t) 


aF(t/a)e‘^^^^ 

(t-a)^/^(t+a) 


(5.3) 


and cr 


o 

yl 


(x,0) 


is a function which is bounded at the end points. 


The behavior of the Cauchy Integral in (5.2) near the end 
points X = ±a will be determined by following the method given by 
Muskhelishvili [8] . 


Hence, defining a sectionally holomorphic function 

a 

/ 

-a 

from (5.3) and (5.4) we can write [8] 


' TT t “ Z 


(5.4) 


(z) = 


/K72 /a72 F(l) ^ 


(z+a) 


1/2 


(z-a) 


1/2 


(5.5) 


Here, is bounded everywhere except the end points where 


$ (z) I < 


I 


< 1/2 


z±a 


C and being real constants. 


Hence from (5.1), (5.2), (5.4) and (5.5) 


K - - 


- V''yx> 


F(l) 


(5.6) 
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which is also the same at the tip x = —a. (5.6) can be derived 

also by using a different approach (see Appendix D) . 

F(l) can be obtained in terms of F(t^), j= l,...,n either by 
an extrapolation or even better by using the formulas given in [11] . 


5.2 For the Material of Type II 

In this case cF^^(XfO) can similarly be found by replacing ^^20 
in (5.2) by [see (4.26)]. Hence, 


K = - 


2A E /a 
21 Y 

IT- V V — r 

xy yx 


F(l) 


(5.7) 


6. DISCUSSION 

From the point of view of the analysis, the important problem 
is the one with nonzero crack surface tractions and vanishing 
loads at infinity. The solution to the actual problem, with the 
stress-free crack surfaces and the loads applied at infinity, can 
be obtained by a simple superposition. Hence, in the following, 
we will assume 

ayi(x,0) = --p(x) ~ -p^ = constant (6.1) 

and plot K/p^/a vs. a/h for various material combinations. 

Two different materials and their elastic constants will be 
given below as examples of materials of Type I . 
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Boron Epoxy 


E = 3.5x10^ 
X 

psi. 

E = 3 

y 

.24x10^ psi 


V =0.23, 

yx 


G 

xy 

1.23X10^ psi 


and for the plane 

strain 



• 

E = 3.5x10^ 
z 

psi. 

V = 

zx 

V = 0.23 

zy 

(6,2) 

Glass-Fiber (20% volume fraction) 



E = 6.6x10^ 
X 

psi. 

E = 2 

y 

.52x10^ psi 


V =0.32, 

yx 


G = 

xy 

2.9x10^ psi 


and for the plane 

strain 




E = 6.6x10^ 
z 

psi. 

V = 

zx 

V = 0.32 

zy 

(6.3) 

In Figure 1, the 

results 

in the case of generalized 

plane 


stress are plotted for two material combinations consisting of 
boron-epoxy and glass-fiber. As may be expected, the lower curve 
illustrates the stiffening effect of a boron-epoxy matrix, result- 
ing in possible crack arrest. Both curves tend to unity as a/h^O 
^bich corresponds to the solution of an infinite orthotropic plate 
with a crack. At the other end, as a/h-^ 1, the effect of the 
material dissimilarity becomes more apparent. If the crack propa- 
gates and reaches the interfaces, then the resulting problem, for 
which a= h, requires a separate treatment. 

Figure 3, like Figure 2, illustrates the variation of K/p^/a 
vs. a/h in the case of plane strain. 

Similar curves are given in Figures 4, 5 and 6 for the 
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materials of Type II. For this purpose the following materials 
are chosen: 


For Figures 4 and 5 : 

(I) : E = 2x10^ psi, E = 2.52x10^ psi 

X y 

V =0.20, = 4x10^ psi 

yx xy 

and for the plane strain 

E = 2x10^ psi, V = V = 0.20 

2 ^ zx zy 

(II) : E = 2.40x10^ psi, E = 8x10^ psi 

X y 

V =0.25, = 3 x 10 ^ psi 

xy xy 

and for the plane strain 

= 2.40x10^ psi, V = V = 0.25 (6.4) 


For Figure 6 : 


(III) : E^ = 2.40x10 psi. 


V “ 0.25, 
yx 


(IV) : E = 1.5x10 psi. 


V = 0.20, 

'^yx 


Ey = 8x10^ psi 

G = 3x10^ psi 
xy ^ 

7 

E =4x10 psx 

y 

G = 4x10^ psi 
xy 


(6.5) 


The discussion of these curves can be made similarly. 
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APPENDIX A 


It can be shown that and S 2 as defined by (3.4) are either 
real or complex conjugates. Hence, two possibilities exist: 


1 ~ Material of type I: 

Sf = ^2 ~ ^3 ^2 “ ^4 ” ^ 

from which it follows that 
^6^ ^7^ ^8^ ^ 9 ' ^10 real 

2 - Material of type II: 

Sj^ = + i(j)2 

S 2 = - i»2 = ii , m^ = 0)3, o)2=-M4 

from which we obtain 

6g = = 6g, ’ 05 = pure imaginary 


(A-1) 


(A-2) 


(A- 3) 


In the case of generalized plane stress from (2.20) and (2.17) 
we have 

A E E 

= t; > - z r , 0n = ^ 


- '’xyV> ' ' " '^xy<l " ^xy V) 


V E 


^3 ^ G (1 - V V ) 

xy' xy yx' 

Hence from (3.4) and (A- 4) 

E 


^4 ^\x “ G ' 

^ xy 


E 

8 r- ~ ' ^ 

^5 E 

X 


(A-4) 


(A-5) 


Therefore the equation to determine and S 2 takes the form [see 
(3.4) ] 
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0 


(A-6) 


^ EE E 

X xy ^ X 


or by change of variables 

i= iu 


we have 


(A-7) 


U 


4 


+ 






0 


(A- 8) 


It has been proved by Lekhnitskii [6] and [12] that equation 
(A-8) has no real roots. The possible roots are 


a) 

= *1^ ' 

U 2 = ' 

Vl , 


or 


b) 

Ml = 63 + 64!, 

U2 = -«3 + «4i . 

^1' 

^2 


(A- 9 ) 


Hence from (A-7) and (A-9) we conclude that the roots in terms 
of s are 


a) 

®1' ®2' 

”®2 

all being real, or 


b) 

S^/ — ^1*^ 

-^1 

s^ = complex 

(A-10) 

I the propositions in 

(A-1) 

and (A- 2) are proved. 


In 

the case of plane 

strain 

the same results are 

obtained [6] , 


[ 12 ]. 
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APPENDIX B 


B, 1 Material Type I ; 


Dimensionless bielastic constants y^z 


Y, = 1 + , 


Y2 = 1 + 


1^3 = “l + V^7 ' 


74 = “3 + V^8 




01-,$. 


Yc 

J. Zf 

= V + , 

Yc = 

V + 


' 5 

/e; 

' 6 

xy 


^7 

= V oIt + 3- , 
xy 1 ' 

= 

V 0)o 
xy 3 

+ 9s 

Y 9 

= -1 + f 

^10 

= -1 + 

0)3$g 

* 

* * 

* 


* * 

79 

= -1 + V 7 ' 

7l0 

= -1 + 

“3®8 


0 ). 


01 . 


Y 


11 


* 

15 


17 


J_ 

+ assign ( 01 ^) , 

^12 

’ ^ 

- + 

3^Qsign(oi3) 

= - hll 

- Vy^B^sign(oij^) 

* 

' *^14 

- I“*3l 

- 


= - 41 

★ I * * 

w^l - B^signCoi^) 

* 

' 7i6 = 

* 1 

- V 

xy ' 

* 1 

“ 3 I 

ic it 

- 3gSign{o)3) 

* 

_ 1 - 

V V 

xy yx 

* 

G 

^ xy 




^X 1 ~ 

* * ' '18 

V V 
w w 

G 

xy 





11 


y^^ = sign(o>j^)Bg - sign{ai3)3^Q 


Y 


Ytr - Y 


11 


6 Y 


Y 


12 


20 


2y 


19 
* * 


* * 


Y21 = sign(oi^)32 " signCo)^)^^ 


40 



Y 


* ^8 


22 = 1 + sign (0)3) ^ 


Y 


^23 “ ^3 


17 


^21 


^7 ^‘^14 ^13^ 


Y 


24 


* * * * ^17^13 

^22^17^14 - sign (< 03 ) 63 


^25 = - + 


17 


^o(Yiii Yio) 


4 " ^8'^14 '13 


* * 


26 


8^<5 

- Yq - -7*— Sign ((03)63 + Y18Y10Y22 


21 


27 


* * 
Yiq ” Yq 


^28 = - Y 


~'^18''''9 
10 ■ Y5i 


* * 


sign(M3)8g + Yj_gYj_o1^22 


Y 


* * 
^10 ■ ^9 

Y.o3o ( — ) 


29 '18^8 


^1 


Y 


Y 


17 


* * 


30 


^21 


^^^13 “ Yi4)3gSign((03) - Y 14 Y 21 I = ”Y 


24 


31 


117 * * 

^ (Yi4 - Y 13 ) 


32 Y 


IP * ** * 

5“ [(Yq - Y^q) 38sign(o)3) - Y;^0^21^ 
2 1 


Y 


Y 


29 


33 6 


^34 = ^30 + YgiBgSignC.^^) 


Y,/^ 


0 ). 


Y 


Y 


11 


35 Y 


12 


tYso + Y3iBioSign((03) - j— ] 


41 



^36 ^32 




(0. 


Yn/?7 


37 

= t^32 

^12 

+ Y33BioSign( 

0 ) 3 ) 

+ 

J. ^ ^ 

l“ 3 l 

] 

38 

^ ^36^25 

^34*^29 

r 

*^39 

= 

*^36^24 

^34^28 

40 

^34‘^27 

"^36^23 

/ 

^41 

= 

^34^26 

^36^24 

42 

”^35^29 

*^37^25 

r 

^43 

5= 

^35^28 

^37^24 

44 

"^37^23 

*^35^27 

t 

^45 

= 

^37'^25 

‘ ^35^26 


B« 2 Material Type II : 


Dimensionless bielastic constants 


= 3^ + i3^ / 


^9 ^9 ^^9 


* * » * « 

3^ = 3^ + i^y , 


* * • * « 

^9 ~ ^9 ^^9 


*1 = “l + V®7 ' 


*2 = “2 + V ®7 


= 


^ ^ ^ M ^ II t 

+ V _^3y f = 1 + {“0)23g + 


yx‘ 




A, = 


> yx 


(o)2sign(o)^) 3g + 


^5 = 
* 

A^ = 


^8 = 

* 

An = 


>- V ^ + 
1 xy 

* * 

3 t 

^6 

2 xy 

I^V + 

3^ , 

A-, = 

V + 

2 xy 


7 

xy 

‘4 

; # 
yx 


Ag = 

-1 - , 


* * »* * * t 

-1 - ^ 2^1 ^ 1^7 


( + ^1^9^ 

/35 


‘2^7 ^1^7 


42 





^10 ' ^ 2^7 ' 


11 




- 3 , 


w. 


12 




+ 0gSign(o)j^) 


13 
* 

14 
* 


= 0) 


lir ^ ^ f ^ 

* ''yx®7 


= v^ylo.;! + e; sign«.') 


* * t 


* . * ” 


^15 “ ^2^7 


‘16 = - egsigtKojj^) - 69 


M ^X2 * 

= - 3^ + ^ — 3gSign(w^) 


12 

18 " ^4 “ ^3 ' 


19 


= ~ A 


11 


12 

11 


‘20 


E (1 - V V ) 
X xy yx 

E 


X (1 - V V ) 
xy yx 


21 


‘2o«; 


‘22 


^20^13 ' 


23 


- xy 

^xy 


‘24 


^23^9 ' 


A 


25 


^23^15 


‘26 


“12 

^7 ^8 ' 


1 ^26 


27 


2 A 


16 


28 


= CO. 


“12 

An ^2 ^ 

11 


29 


= ll (|o)^| 


+ (0 




2 A 


ii 

11 


‘12 


30 


All 


31 


• 1 ^30 

■ 2 Ai, 
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29 


w, 


32 


^30“l "l ' 


^4 


1 

^^33 

2“2*16> 

' ^35 

A 



^34 


A s= 

^36 


^1 

*35“l 

r 

37 

A 


1 ' 



A = 

^38 


2 

Ai 6 

r 

^39 



o 

^^21^7 

*ii 


^40 


A 39 

- *1^7 

) . = 



9 

t 

* • 

A 

^42 

“ 

^39 

^^22^7 

“ A^3^ 

sign ((D^) ] 



9 

II 

* » 

* 

A > 
43 

““ 

39 

(*22®7 

- *2^7 

sign (w^) ] 

A 


1 

1*40 (* 

24 "^9^ 

— A-» A * « 

^44 


^^10 

25 42 

^45 


1 

■ I*4i(* 

24 "^9^ 

— A _A.^ 


2*10 

25 43 

*46 


Ag 


A = 

1 r 


^^10 

t 

^47 

2 A A ^ 

^^10 39 


^33 ~ “2^16 “ l“ll *17 

_ *18*^ . 

2(A), A-, ^ 

1 16 


*'» 


* » 


•y H * »» 


‘24” 9^^21 25 22 


48 


2Art t- + A256*'sign(«)*)] 

10 39 


*11 


‘49 


‘21 


2A ' 


'50 2A 
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51 


^^ 31 ^ 32^46 Aj_q 0 )j_ ^47 


‘52 


= A^^A^ , 


^53 “ ^51^6 


‘54 


” ^31^32 ■*’ ^50 ' 


^55 “ ^54^5 
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56 




^38 

'57 + A3^A4q) + ~ ^36^47 


'58 


^ 57^5 ' 


^59 - ^ 57^6 


‘60 “ ^31 ^36^50 ~ ^^37^49 


‘61 


^ 60‘^5 ' 


^62 “ ^ 60^6 


63 


OA A A A 

^^ 31 ^ 32^46 “ ^47 " A ^^ 


'64 


^ 63^6 ' 


A^c = “ Ac-^Ap. 

65 63 5 


'66 


^ 31^32 ^50 


‘67 


" * 66*6 ' 


Aco - A-^Ap. 
68 66 5 


69 

^^^ 31^46 

70 "" 

^ 69^6 ' 

72 "" 

- ^31 - 

73 "" 

"' ^ 12^6 


A^q 36 47 


^71 ““ ■“ ^ 69^5 


^ 49^37 


^74 ” ^ 72^5 
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APPENDIX C 


Calculation of certain integrals « 


From [9] we have 


/ e ^^sin bt dt = 


2^,2' 
a + b 


a > 0 


/ e ^^cos btdt = 


2 k 2 ' 

a + b 


a > 0 


(C-1) 


/ 

o 


+ (Y-x)^ 8^ + (Y+x)^J 


sin(ax)dx = J e ^®sin(aY) > 


o [e^ + (Y-x)^ ' 8^ + (Y+x)^J 


a> 0, Re(B) > 0, (Y+i3) is not real 
1 


cos(ax)dx = ^ e ^^cos(ay) r 


a > 0 , I Imy | < Re ( 3 ) 

and by properly differentiating with respect to a 


X 


3^ + (Y“X) ^ 3^ + (y+x)^ 


cos (ax) dx 



TT 


“ F ' 

CO 

_ 

J 

o 

x 

.8 ^ + (y 


sin (ax) dx 


= j [Ysin (ay) + 3cos(ay)] 


Special Cases , 


(C-2) 


DO 

/ ^ ^ — o cos(ax)dx = e ^ 


o2 ^ 2 

o 3 + X 


23 


/ 2 ^ — 2 “ 5 ^ 

o + x'^ 
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/ ^ sin(ax)dx = (1 - e 

o x(3 +x^) 23 


a + b/^ 

+ (d^+ ft)^ 


b(c + d^) ^ 2bd b 

^2^2 t, + a - 2 2 

f t ^ t 

cC^ + (d^ +ft)^ C 


(C-3) 
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APPENDIX D 


Derivation of the Stress Intensity Factor Expression: 


From (5-2) and (5.3) we have 

1 




(ax,0) 


= w(l - v,;v -)- _[ - -yl 

where x = ax / t = aT , 


(D-1) 


Expanding F(t) into a series of Chebyshev polynomials 

00 

F(t) = I C.T.(t) (D-2) 

j=0 ^ ^ 

and substituting into (D-1), with the use of the formula [13] 

f 1x1 < 1 

1x1 > 1 (D“3) 




= < 


-1 /rr^ ^ 

for Ixl > If obtain 


tt[/y^-1 ~ X]^ 

I (-1)3+1 ' 


cJyl(aXfO) = 


^~ ^20^y 


I C. 1 + (aXfO) , 

(1 - v^yVy^)/^ j=0 3 (_1)3 


yi 

Ixl > 1 (D-4) 


Hence 


K =. lim /2 (x-a) cr n (x,0) = 
x^a ^ 


^1 " '^xy'*yx* j=0 3 


(D-5) 


But from (D-2) 


oo 

F(l) = I C.T.(l) = I C. 
j=0 j=o 


(D-6) 


Hence, we finally obtain 


K = - 


2y2oK'^ 

y - - F(i) 


” '*xy'^yx^ 


(D-7) 
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Fxgur© 2« K/pi^s vs* a/h v 

Material of Type I, Generalized Plane Stress 




Figure 3. K/p/a vs. aA 

Material of Type I 
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Plane Strain 






